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Abstract. We investigate higher-order Voronoi diagrams in the city metric. This metric is 
induced by quickest paths in the L\ metric in the presence of an accelerating transportation 
network of axis-parallel line segments. For the structural complexity of fc th -order city Voronoi 
diagrams of n point sites, we show an upper bound of 0(k(n — k) + kc) and a lower bound 
of 0(n + kc), where c is the complexity of the transportation network. This is quite different 
from the bound 0(k(n — k)) in the Euclidean metric [11]. For the special case where k — n— 1 
the complexity in the Euclidean metric is 0(n), while that in the city metric is 0(nc). 
Furthermore, we develop an 0(k 2 (n + c) logn)-time iterative algorithm to compute the k th - 
order city Voronoi diagram and an O (nc log 2 (n + c) logn)-time divide-and-conquer algorithm 
to compute the farthest-site city Voronoi diagram. 

1 Introduction 

In many modern cities, e.g., Manhattan, the layout of the road network resembles a grid. 
Most roads are either horizontal or vertical, and thus pedestrians can move only either 
horizontally or vertically. Large, modern cities also have a public transportation network 
(e.g., bus and rail systems) to ensure easy and fast travel between two places. Traveling in 
such cities can be modeled well by the city metric. This metric is induced by quickest paths 
in the L\ metric in the presence of an accelerating transportation network. We assume that 
the traveling speed on the transportation network is a given parameter v > 1. The speed 
while traveling off the network is 1. Further, we assume that the transportation network 
can be accessed at any point. Then the distance between two points is the minimum time 
required to travel between them. 

For a given set S of n point sites (i.e., a set of n coordinates) and a transportation 
network in the plane, the k^-order city Voronoi diagram Vk(S) partitions the plane into 
Voronoi regions such that all points in a Voronoi region share the same k nearest sites with 
respect to the city metric. 

The A: th -order city Voronoi diagram can be used to resolve the following situation: a 
pedestrian wants to know the k nearest facilities (e.g., k stores, or k hospitals) such that he 
can make a well-informed decision as to which facility to go to. For this kind of scenario, 
the /c th -order city Voronoi diagram provides a way to determine the k nearest facilities, by 
modeling the facilities as point sites. 

The nearest-site (first-order) city Voronoi diagram has already been well-studied |1|4|6|I10| . 
Its structural complexity (the size) has been proved to be 0(n + c) jl], where c is the 
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complexity of the transportation network. Such a Voronoi diagram can be constructed in 
0((n + c) log(ra + c)) time [6]. However, to the best of our knowledge there is no existing 
work regarding /c th -order or farthest-site (i.e., (n — l) st -order) city Voronoi diagrams. 

Contrary to fc th -order city Voronoi diagrams, /c th -order Euclidean Voronoi diagrams 
have been studied extensively for over thirty years. Their structural complexity has been 
shown to be 0{k{n — k)) [11 . They can be computed by an iterative construction method 
in 0{k 2 n\ogn) time [11] or by a different approach based on geometric duality and ar- 
rangements in 0(n 2 + k(n — k) log 2 n) time [7]. Additionally, there are several randomized 
algorithms |2|13j and on-line algorithms [3)3]. 

One of the most significant differences between the Euclidean metric and the city metric 
that influences the computation and complexity of Voronoi diagrams is the complexity of 
a bisector between two points. In the Euclidean or the L\ metric such a bisector has 
constant complexity, while in the city metric the complexity may be Q(c) [T] and can even 
be a closed curve. Since the properties of a bisector between two points significantly affect 
the properties of Voronoi diagrams, a /c th -order city Voronoi diagram can be very different 
from a Euclidean one. First, this property makes it non-trivial to apply existing approaches 
for constructing Euclidean Voronoi diagrams to the city Voronoi diagrams. Secondly, this 
property also indicates that the complexity of /c th -order Voronoi diagrams may depend 
significantly on the complexity of the transportation network. 

In this paper, we derive bounds for the structural complexity of the fc th -order Voronoi 
diagram and develop algorithms for computing the A: th -order city Voronoi diagram. The 
remainder of this paper is organized as follows. In Section [2j we introduce two important 
concepts, wavefront propagation pQ and shortest path maps [B], which are essential for the 
proofs in the subsequent sections. In Section [3j we adopt the wavefront propagation to 
introduce a novel interpretation of the iterative construction method of Lee [llj, and use 
this interpretation to derive an upper bound of 0(k(n—k)+kc) for the structural complexity 
of A; th -order city Voronoi diagrams, where c is the complexity of the transportation network. 
Then, we construct a worst-case example to obtain a lower bound of i?(n + kc). Finally, 
we extend the insights of Section [3] to develop an iterative algorithm to compute k th - 
order city Voronoi diagrams in 0(k 2 (n + c) log(n + c)) time (see Section [4]). Moreover, 
we give a divide-and-conquer approach to compute farthest-site city Voronoi diagrams in 
0(nc log n log 2 (n + c)) time. We conclude the paper in Section pi 



For an overview of our contribution and a comparison between Euclidean and city 
metric see Table [U 

2 Preliminaries 

In this section we introduce the notation used throughout this paper for fc th -order city 
Voronoi diagrams. Then, we introduce two well-established concepts in the context of 
Voronoi diagrams, which are important for the proofs in the subsequent sections. 

A transportation network is a planar straight-line graph C = (Vc,Ec) with isothetic 
edges only, i.e., edges that are either horizontal or vertical, and all transportation edge have 
identical speed v > 1. We define c := \ Vc\, and since the degree of a vertex in Vq is at most 
four, \Ec\ is 0(c). We denote the distance of two points in the L\ metric by d\ and in the 
city metric by dc- Similarly, we denote the bisector between two points by B\ and Be for 
the L\ and city metric, respectively. Additionally, for the city metric we define the distance 
between a point p £ M 2 and a set of points H C M 2 to be dc(p,H) = max 5£ ^ dc(p, H). 
This allows us to define the bisector Bc{H\, H2) = {r £ I 2 dc{r,H\) = dc(r,H2)} 
between two sets of points Hi and H2. 

By Vk(H,S) we denote a Voronoi region of Vk(S) associated with a /c-element subset 
H C S. The common boundary between two adjacent Voronoi regions Vk(Hi,S) and 
Vk(H2, S) is called a Voronoi edge. This Voronoi edge is a part of Be (Hi, H2) = Bc{p, q) 
where H\ \ H2 = {p} and H2 \ Hi = {q} The common intersection among more than 
two Voronoi regions is called a Voronoi vertex. Without loss of generality, we assume that 
no point in the plane is equidistant from four sites in S with respect to the city metric, 
ensuring that the degree of a Voronoi vertex is exactly three. 

Wavefront Propagation. The wavefront propagation is a well-established model to de- 
fine Voronoi diagrams P. In Section g we will use this concept to interpret the formation 
of Vk(S) and analyze its structural complexity. 

For a fixed site p 6 S, let W p (x) = {q \ q G M. ,dc(p, q) = x}. This means that for 
a fixed x £ Wq the wavefront W p (x) is the circle centered at p with radius x. We call p 
the source of W p (x). Note that we can view W p (x) as the wavefront at time x of the wave 
that originated in p at time 0. We refer to such a wavefront as W p if the value of x is 
unimportant. 

Initially, the wavefront W p is a diamond. When it touches a part of the transporta- 
tion network for the first time it changes its propagation speed and, hence its shape; see 
Fig. [T| Certain points on the transportation network play an important role to determine 
the structural complexity of fe th -order city Voronoi diagrams. Thus, we introduce the fol- 
lowing definitions. For a point v £ M 2 , let P(v) denote the isothetic projection of v onto 
the transportation network, i.e., we shoot an isothetic half-ray starting at v in each of the 
four directions and for each half-ray we add its first intersection with an edge of the trans- 
portation network to P(v). It is easy to see that there are at most four such intersections. 




Fig. 1. Wavefront Propagation. 



Fig. 2. 1-needle, 2- needle and 3-needle. 
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Fig. 3. (a) A 2-needle is two 1-needles. (b) Wavefront propagation of a 1-needle. (c) Bi(rj p (pi), iq q (q\)). 

For a set I C I 2 , we denote the isothetic projection of the set X as V{X) = y} v&x P( v )- 
For a site p £ S, we call the set A(p) = P(p) U Vc U V{Vc) U {p} activation points (we 
added {p} to the list for ease of argumentation in some of our proofs). 

As shown by Aichholzer et al. pQ, the wavefront W p changes its propagation speed only 
if it hits a vertex in A(p). Since the shape of W p can become very complex after it hits 
multiple activation points, we make the following simplification for the remainder of this 
paper: if a wavefront W p touches a point q E A(p) we do not change the propagation speed 
of W p . Instead, we start a new wavefront at q, which, in turn, starts new wavefronts at 
points in A{p) if it reaches them earlier than any other wavefront. Hereafter, the start of 
the propagation of a new wavefront is called an activation event, or we say a wavefront 
is activated. The shape of such a new wavefront depends on the position of q on the 
transportation network. It can be categorized into one of three different shapes: 1-needle, 
2-needle, and 3-needle [Ij (see Fig. [2]). To simplify things, we treat a 2-needle (3-needle) as 
two (three) 1-needles (see Fig. [3j[a)). 

When a 1-needle reaches the end of the corresponding network segment, as shown in 
Fig. [3^b), its shape changes (permanently) [I]. In order to interpret the propagation of a 
1-needle, Bae et al. [6] introduced a structure called needle. A needle r) p (q, q') is a network 
segment qq' with weight dc(p,q), where p £ S and q,q' £ A(p). Propagating a wavefront 
from r]p(q,q') is equivalent to propagating a 1-needle from q on the network segment qq' 



at time dc{p,q)- If q' is obvious or unimportant we may refer to r] p (q,q') as r] p (q). Bae et 
al. also defined the L\ distance di(rj p (q, q'),r) between a needle rj p (q,q') and a point r as 
dc(p, q) plus the length of a quickest path from q to r accelerated by qq'. Thus, the bisector 
Bi(r]p(pi), r)q(qi)) between two needles r) p (pi) and r) q (qi) is well defined (see Fig. [3^c)). 

Shortest Path Map. We use the wavefront model to define shortest path maps [6], 



and use this concept to explain the formation of mixed vertices in Section 3.1, which are 
important for deriving the structural complexity of the fc th -order city Voronoi diagram. 

For a site p G S its shortest path map SVM. P is a planar subdivision that can be 
obtained as follows: start by propagating a wavefront from the site p. When a point q £ A(p) 
is touched for the first time by a wavefront, propagate an additional wavefront from rj p (q). 
Eventually, each point r £ M. 2 is touched for the first time by a wavefront propagated from 
a needle r] p (q), where q G A(p) and di(r,r) p (q)) = min q i e ^ p -jdi(r,r] p (q')), and q is called 
the predecessor of r. This induces SVM. P . In detail, SVM. P partitions the plane into at 
most |^4(p)| = 0(c) regions SVM. p (q) such that all points r G SVM. p {q) share the same 
predecessor q and q is on a quickest path from p to r, i.e., dc(p, r) = dc(p, q) + dc(9> r ) = 
di(r,r] p (q)). As proved in [6], the common edge between SVM. p {q) and SVAi p (q') where 
q, q' G ^4(p) belongs to the bisector Bi{rj p {q),rj p {q')). Fig. [4] illustrates an example of the 
function of shortest path maps where the two Voronoi regions of Vi({p, q}) are partitioned 
by SVM. P and SVA4 q , respectively. 

3 Complexity 

In this section we derive an upper and a lower bound of the structural complexity of the 



/c th -order city Voronoi diagram Vk(S). In Section 3.1 we first introduce a special degree-2 
vertex on a Voronoi edge called mixed vertex which is similar to the mixed Voronoi vertices 
of Cheong et al. [8] for farthest-polygon Voronoi diagrams. Then we derive an upper bound 
of the structural complexity of (S) in terms of the number of mixed vertices and Voronoi 



regions. In Section 3.2, we adopt the wavefront concept to introduce a new interpretation 



for the iterative construction of Vk(S) by Lee [TT]. This yields an upper bound for the 



structural complexity of Vk(S). In Section 3.3 we construct a worst-case example to obtain 



a lower bound for the structural complexity of Vk(S). 



3.1 Mixed Vertices 

Definition 1 (Mixed Vertex) For two sites p,q G S and a Voronoi edge e which is part 
of Bc(p, q), a point r on e is a mixed vertex if there are pi,P2 G A(p) and q\ G A(q) such 
that r G STM p (pi) n STM p {p 2 ) n STM q (qi). 

For instance, Fig. [4] shows a first-order city Voronoi diagram V\({p, q}), where the 
mixed vertices are marked with a square and denoted by mi, . . . , 711,4. The vertex 7712 is a 




Fig. 4. Bc{p, q) (solid thin edge), where mi, . . . , TO4 are mixed vertices. 

mixed vertex because it is in SVMp{pi)r\SVM p (p 2 )nSVMq{qi). Definition [l] yields the 
following. 

Lemma 1 If a Voronoi edge e contains m > mixed vertices, its complexity is 0(m + 1). 

Proof. Suppose e is part of a bisector Bc(p,q), where p,q E S. Consider two consecutive 
mixed vertices mi and vri2 on Bq(p, q), where m\ E SVM. p (p\) n SVM. q (q' 1 ) n 5"PA^ g ((/i) 
and m 2 E SVM p {pi) n SVM p {p2) n SP-M^i) (see Fig. [3]). Consider each point u 
on Bc(p,q) between mi and 7712- Since v belongs to SVM p (pi) n SVAi g (qi), we have 
d c {v,p) = d c (v,pi) +d c (pi,p) = di{v,ri p (pi)) and d c (v,q) = d c {v,qi) + d c (qi,q) = 
di(v,r) q (qi)). Together with d c (v,p) = d c (v,q), v belongs to B^rjpipi) , 77 9 (<?i)) (recall 
Figj^c)). As a result, if a Voronoi edge e contains m mixed Voronoi vertices, e consists of 
m + 1 parts, each of which belongs to a bisector between two needles. Since the complexity 
of an L\ bisector between two needles is 0(1) [3], the complexity of e is 0(m + 1). Note 
that the complexity of a bisector between two points in the city metric is 0(c), while the 
complexity of an L\ bisector between two needles is O(l). □ 

Lemma 2 An upper bound for the structural complexity of a k til - order city Voronoi dia- 
gram Vk(S) is 0(M + k(n — k)), where M is the total number of mixed vertices. 

Proof. Lee [11] proved that the number of Voronoi regions in the /c th -order Voronoi diagram 
is 0(k(n — k)) in any distance metric satisfying the triangle inequality, and so is the number 
of Voronoi edges. By Lemma [TJ if a Voronoi edge e contains m e mixed Voronoi vertices, its 
complexity is 0(m e + l). Suppose a city Voronoi diagram Vk(S) contains a set E of Voronoi 
edges, and each edge e E E contains m e mixed Voronoi vertices. Then, the complexity of 
all edges, i.e., the structural complexity of Vfc(S), is X^ee-E 0(m e + 1) = 0(M + \E\). Since 
\E\ = 0(k{n - k)), it follows that 0{M + \E\) = 0{M + k(n — k)). □ 



Fig. 5. Vi(H,S)nVi(Q) = Vi(H,S)nV 2 (S) where Fig. 6. V 2 (H,S) where H = {pi,p 2 }, Q = 

H = {p}, Q = Ui<i< 6 {%}» and 5 = HUQ. { qi ,q 2 ,q a ,q 4 }, and S = U Q. 



For the proof in Section 3.2 we further categorize the mixed vertices. Let m be a mixed 
vertex on the Voronoi edge between Vk(Hi,S) and Vk(H2, S), where H\ \ H2 = {p} and 
H2 \ H± = {q}. We call m an interior mixed vertex of Vk(Hi,S) if m 6 SVM p (pi) D 
SVM. p {p2) H SVAi q (qi), for some pi,P2 £ -^(p) and (/1 E A(q); otherwise, we call m an 
exterior mixed vertex of Vk(Hi, S). For example, in Figj4]the vertices ni2 and m4 both are 
interior mixed vertices of Vi({p}, {p, q}) and exterior mixed vertices of Vi ({</}, {p, q}). 



3.2 Upper Bound 

Throughout this subsection, we consider a Voronoi region Vj(H, S) of a j th -order Voronoi 
diagram Vj(S), where H C S and \H\ = j. Let Vj(H, S) have adjacent Voronoi regions 
Vj(Hi, S) for 1 < i < hu. Note that the subsets Hi and H differ in exactly one element 
In the following let Hi\H = {%}, Q = {<?i, . . . , qh H }, and ^ = \Q\. 

Lee [ll] proved that in any distance metric satisfying the triangle inequality, Vj (H, S) n 
Vi(Q) = Vj(H, S)nVj + i(S), and thus computing V\{Q) for all the Voronoi regions Vj(H, S) 
of Vj(S) yields Vj + \(S), leading to an iterative construction for Vf.(S) for any k < n. Fig. [5] 
illustrates this iteration technique for the Euclidean metric: solid segments form V\(H, S) 
and dashed segments form V\(Q). Since the gray region is part of Vi({p}, S) and also part 
of V\({qi},Q), all points in the gray region share the same two nearest sites p and qi, 
implying that the gray region is part of V2({p, qi},S). 

We adopt wavefront propagation to interpret this iterative construction in a new way, 
which will lead to the main proof of this section. Let us imagine that a wavefront is 
propagated from each site q £ Q into the Voronoi region Vj(H, S). If a point r £ M. 2 is first 
touched by the wavefront that propagated from q, r belongs to Vj + \{H U {q}, S). 

Note that when j > 2, \Q\ is not necessarily the number of adjacent regions, i.e., In < 
htf. Fig. [6] illustrates an example for the Euclidean metric: V^-ff, S) has 6 adjacent Voronoi 



regions but \Q\ = £h is only 4. This is because for a site q £ Q, B2({q}, H) n Vj(H, S) may 
consist of more than one Voronoi edge, where B2({q}, H) is a Euclidean bisector between 
{q} and iif (similar to Bc(Hi, H2) defined in Section[2]). For instance, as shown in Fig. |6j 
e qi = £?2({<7i}i H) n V2(H, S) consists of two Voronoi edges ei and e2 

Now we transfer our new interpretation to the city metric. Let e q be Bc({q}, H) n 
Vj(H,S) for some site g G Q. If e g contains m g exterior mixed vertices with respect to 
Vj(H,S), e q intersects m q + 1 regions in SVM. q . We denote these regions by SVAi q {v z ) 
for 1 < z < m g + 1. Note that all u 2 must be in A(q). Then, instead of propagating a 
single wavefront from q into Vj(H,S) (as in the Euclidean metric), we propagate m q + 1 
wavefronts, namely one from each r] q (v z ) into Vj(H,S). 

As a result, if Vj(H, S) contains ran exterior mixed vertices, ran + £h wavefronts will 
be propagated into Vj(H, S). During the process, when a point r £ Vj(H, S) is first touched 
by a wavefront propagated from n q (v),q £ Q and v £ A(q), we propagate a new wavefront 
from rj q (r), i.e., an activation event occurs, if (i) r £ ViVc) U Vc or (ii) u = q and r £ P(g). 
These two conditions amount to r £ A(q r ) \ {g}, but this classification will help us to derive 
the number of mixed vertices. This is due to the fact that during the k — 1 iterations for 
computing Vj + i(S) from Vj(S) for 1 < j < k — 1, V(Vc) contributes 0(kc) activation 
events, but V(S) only contributes O(n). 

Lemma 3 IfVj(H,S) contains run exterior mixed vertices, then Vj(H, S) P\Vj+i(S) con- 
tains at most run + 2c# + 2an mixed vertices, where ch = \ ('P(Vc) U Vc) H Vj(H, S)\ and 
an is the number of activation events associated with points in V(S). 

Proof. According to the above discussion, we propagate tuh + C-h wavefronts into Vj(H, S). 
All those wavefronts combined generate at most ch new wavefronts from points in 'P(Vc , )n 
Vj{HS), and an new wavefronts from points in V(S)r\Vj(H, S). Note that ch = \(T(Vc)^ 
V c ) n Vj(H,S)\ (condition (i)) but a H < \V(S) n Vj(H,S)\ (condition (ii)) . Let W be 
the set of the run + £h + ch + an wavefronts. For each point r £ Vj(H, S), if r is first 
touched by a wavefront w £ W it is associated with w. This will partition Vj(H, S) into 
m-H + + ch + an regions. We view those regions as a special Voronoi diagram Vi(VF). 
Note that mn + £h of those regions are unbounded. 

Vj(H, S) n Vj+i(S) is a subgraph of Vi(W) since if a point r £ Vj(H, S) is first touched 
by a wavefront in W propagated from rj q (v) , r belongs to Vj + i(HU{q}, S). Without loss of 
generality, we assume every vertex of Vi(W) has degree 3. According to Euler's formula it 
holds that Ny = 2(Nr—1)—Nu, where Ny, Nr and Nu are the numbers of vertices, regions, 
and unbounded regions, respectively. Since V\ (W) contains mu + 0-h unbounded regions 
and mn+^H + CH + ciH bounded regions, Vi(W) contains mn+iii + ^CH -\-2ajj — 2 vertices. 
By [11) . since \Q\ = £h, there are ig — 2 Voronoi vertices in Vj+i(S) n Vj(H, S). Therefore, 
Vj(H, S)nVj + \(S) contains at most (itih+£h + 2ch + 2oh — 2) — (£h — 2) = rriH + 2cH + 2aH 
mixed vertices. □ 



Applying Lemma [3] to each region of Vj(S), yields a recursive formula for the total 
number of mixed vertices rrij+i in Vj+\(S): rrij+i = rrij + 0(c) + aj (see Lemma [4]). In 
Lemma [5] we show that this formula can be bounded by 0(n + kc) for k iterations of 
this iterative approach. Finally, in Theorem [T] we combine the insights of Lemma [2] and 
Lemma [5] to give an upper bound for the structural complexity of Vk(S). 

Lemma 4 Vj+\(S) contains rrij + 0(c) + 2a j mixed vertices where rrij is the number of 
mixed vertices ofVj(S) and aj is the number of activation events associated with points in 
V(S) during the computation ofVj + \(S) from Vj(S). 

Proof. For a Voronoi region Vj(H, S), let mn be the number of its exterior mixed vertices, 
let ch be \ Vj(H, S) n (P(Vc) U Vq)\ and let an be number of activation events associated 
with vertices in V(S) n Vj(H, S) during the computation. If Vj(H, S) is empty, m# = ch = 
an = 0. By Lemma|3j Vj(H, S)C\Vj+i(S) contains at most mn + 2ch + 2au mixed Voronoi 
vertices. Therefore, the total number of mixed vertices of Vj+i(S) is bounded by: 

(m H + 2c^ + 2a H ) = rrij + 2\V{V G )\ + 2^. 

H&S,\H\=j 

□ 

Lemma 5 The number of mixed vertices of Vk(S) is 0(n + kc). 

Proof. Let rrij be the total number of mixed Voronoi vertices of Vj(S) and let aj be the 
number of activation events associated with vertices in 'P(S) during the computation of 
Vj+i(S) from Vj(S), described by our algorithm. Then, by Lemma [4] the following holds: 

mk = mfc_i + 0(c) + au-i = ■ ■ ■ = m\ + 0(kc) + 2 ^ aj. 

i<i<fe-i 

Now, we show an upper bound for the complexity of J2i<j<k~i a r F° r a vertex v S P(q) 
where q G S, let the j' th iteration be the first time when v is activated by a wavefront 
propagated from q, i.e. rj q (v) will propagate a wavefront, and let v belong to Vj(H, S). Due 
to this and since the points in H are the j nearest sites of v, q is the (j + l) st nearest site 
of v. Therefore, for j' > j, if v E Vji(H' , S), q G H' , implying that v will not be activated 
by a wavefront propagated from q again after the j th iteration. In other words, v causes at 
most one activation event due to the wavefront propagation of q during the k — 1 iterations, 
and the V(S) causes 0(n) activation events, i.e., J2i<j<k-i a j = 0(n). Furthermore, mi 
has been proved to be 0(n + c) |l|6lll0j . Therefore, m/% = 0(n + kc). □ 



Theorem 1 The structural complexity ofVk(S) is 0(k(n — k) + kc). 
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Fig. 7. This worst-case example (here with k — 3, n = 12, c = 18) leads to a lower bound of Q(n + he). 
The bold solid segments depict the transportation network, and the dashed segments compose Vk(S). The 
right part is also the farthest-site city Voronoi diagram of {s\, S2, S3, S4}, where all points in Region i share 
the same farthest site Si. 

3.3 Lower Bound 

We construct a worst-case example (see Fig. to derive a lower bound for the structural 
complexity of the /c th -order city Voronoi diagram Vk(S). The example consists of a left part 
and a right part which are placed with a sufficiently large distance between them. We place 
one vertical network segment in the left part and build a stairlike transportation network 
in the right part. Then, we place k + 1 sites in the right part and the remaining n — k — 1 
sites in the left part. Since the distance between the left and right part is extremely large, 
the n — k — 1 sites in the left part hardly influence the formation of Vk(S) in the right part. 
Therefore, Vk(S) in the right part forms the farthest-site city Voronoi diagram of the k + 1 
sites, because sharing the same k nearest sites among k + 1 sites is equivalent to sharing 
the same farthest site among the k + 1 sites. 

By construction, as shown in the right part of Fig. [7| all the points in Region i share 
the same farthest site Sj. Since we can set the speed v to be large enough, for each point 
x in Region 2, the shortest path between x and s\ ($2) moves along the transportation 
network counterclockwise, and thus dc(s2,x) > dc(s±,x). The common Voronoi edge be- 
tween Regions i and (i + 1) contains at least (^jp — 1) • 2 + 1 (here: 7) segments since 
the transportation network forms rectangles and each rectangle except the first one 
contains two vertices of the Voronoi edge. Therefore, in the right part, Vk(S) contains at 
least (k — 1)— ?p = Q{kc) segments. Together with the f2(n — k) in the left part, we obtain 
the following lower bound. 

Theorem 2 The structural complexity ofVk(S) is Q{n + kc). 

Proof. We need to distinguish two cases: 

i) n — k — l>k + l=^2k<n — 2: This implies that the left part contains n — k — 2 
segments, and thus, Vk(S) contains n — k — 2 + Q{kc) = Q(n + kc) segments. 



ii) 2k > n — 2: This implies that the left part is empty, and thus, Vk(S) contains 
Q{kc) = Q(k(c + 2)) = Q(n + kc) segments. 

This concludes the proof. □ 



4 Algorithms 

In this section we present an iterative algorithm to compute fc th -order city Voronoi diagrams 
in 0(k 2 (n + c) log(n + c)) time. Its main idea has already been introduced in the complexity 
considerations in Section |3.2| For the special case of the farthest-site Voronoi diagram, 
i.e., the (n — l) st -order Voronoi diagram, this algorithm takes 0(n 2 (n + c) log(n + c)) 
time. However, for the farthest-site city Voronoi diagram we present a divide-and-conquer 
algorithm which requires only 0(nc log 2 (n + c) logn) time. 



4.1 Iterative Algorithm for fc th -Order City Voronoi Diagrams 

We describe an algorithm to compute A; th -order city Voronoi diagrams Vk(S) based on 



the ideas in Section 3.2 and Bae et al.'s j6] 0{(n + c) log(n + c))-time algorithm for the 
first-order city Voronoi diagram V\{S). Bae et al.'s approach views each point site in S as 
a needle with zero-weight and zero-length, and simulates the wavefront propagation from 
those needles to compute V\{S). Since their approach can handle general needles, we adopt 



it to simulate the wavefront propagation of Section 3.2 to compute Vj+i(S) from Vj(S) 



Algorithm. We give the description of our algorithm for a single Voronoi region Vj(H,S). 
All four steps have to be repeated for each Voronoi region of Vj(S). 

Let Vj(H, S) have h adjacent regions Vj(Hi, S) with Hi \ H = {q{\ for 1 < i < h and 
let Q = Ui<i</i#i- Our algorithm computes Vj(H, S) n Vj+i(S) as follows: 

1. Compute a new set N of sites (needles): For 1 < i < h, if the Voronoi edge between 
Vj(Hi,S) and Vj(H,S) intersect mi regions SVM. qi (v z ) in SVM. qi , 1 < z < rm, insert 
every r] qi {v z ) into N. 

2. Construct a new transportation network Ch from C: For each point v € {ViVc) U 
V(Q) U Vc) n Vj(H, S), if u is located on an edge e of C, insert e into Ch- 

3. Perform Bae et al.'s wavefront-based approach to compute Vi(N) under the new trans- 
portation network Ch- The approach can intrinsically handle needles as weighted sites. 

4. Determine Vj(H, S) n V j+1 (S) from Vi(N); Consider each edge e in Vj(H, S) n Vi(N). 
Let e be an edge between Vi{rj p (v p ), N) and Vi(r] q (v q ), N) where p,q £ S, v p G A(p) 
and v q £ A(q). lip^q, then e n Vj(H, S) is part of Vj(H, S) n V j+ i(S). 

Note that Step 2 is used only to reduce the runtime of the algorithm. Lemma [6] shows 
the correctness and the run time of this algorithm for a single Voronoi region. 

Lemma 6 Vj(H, S)D V ? +i(S l ) can be computed in 0((/i + m + c#) log(ra + c)) time, where h 
is the number of Voronoi edges, m is the number of mixed vertices, and ch = {(^(Vc) U 

v c )nVj(H,s)\. 



Proof. We begin by proving correctness. Since Vj(H, S)(lVj+i(S) is exactly Vj(H, S)P\Vi(Q) 
[TT] . it is sufficient to prove that the algorithm correctly computes Vj(H,S) H V\(Q). If 
the algorithm fails to compute Vi(q,Q) n Vj(H, S),q 6 Q, it must fail to propagate a 
wavefront from a needle n q (v), where v belongs to A(q) and SVM. q (v)C\ V\(q, Q)(lVj(H, S) 
is nonempty. We prove that this cannot occur by contradiction. Assume that the algorithm 
does not propagate a wavefront from an rj q (v) for some v in A(q) and SVM. q (v)C\ V±(q, Q)d 
Vj(H,S) is nonempty. However, either v £ Vj(H,S) and SVM q (v) n B c {q,p) D Vj-(iT,S) 
must be nonempty, then, Step 1 will include r] q (v) in iV. Or v £ Vj(H,S), then Step 2 
will include the corresponding network segment in Ch, and thus rj q (v) will be activated to 
propagate a wavefront. Both possibilities contradict the initial assumption. Therefore, the 
algorithm correctly computes Vj(H,S) D Vj+i(5). 

We proceed by giving time complexity considerations. It is clear that \N\ is 0(m + h). 
The run time of step 1 is linear in the complexity of the boundary of Vj (H, S) and thus 
is 0(m + h). Since by definition \{V{V C ) U V c ) n tf/(#,5)| = O(cjy) and \V{Q)\ = 0(h), 
both |Vc H | and \Eq h \ are in 0(c# + h). Since |JV| = 0(m + /i) and -Ec*h = 0(ch + /i), 
Step 3 takes 0((h + m + ch) log(/t + m + ch)) time [6]. Step 4 takes the time linear in the 
complexity of Vi(N) n V}+i(S'). The activation events associated with vertices in V(S) are 
only associated to vertices in V(Q), we know that ^(Q)] = OQi). Therefore, since there 
are 0(m + h), 0(c H ), and 0(h) wavefronts due to N, (V(V C ) U Vc) n ^(F, 5), and P(Q), 
respectively, the complexity of Vi(N) is 0(m + h) + O(ch) + 0(h) = 0(h + m + ch). Since 
m = 0(n + kc) = 0(nc) it holds that 0(\og(h + m + ch)) = 0(log(nc)) = 0(log(n + c) 2 ) = 
0(log(n + c)). We conclude that the total running time is 0((h + m + ch) log(n + c)). □ 

Applying Lemma[6]to each region of Vj(S) combined with Theorem[l]leads to Lemma 
The summation of 0((j(n—j)+jc) log(re+c)) in Lemma[7]for 1 < j < k—1 gives Theorem 

Lemma 7 Vj+x(S) can be computed from Vj(S) in 0((j(n — j) + jc) log(n + c)) time. 

Proof. For a Voronoi region Vj(H, S), let hn be the number of Voronoi edges, m' H be the 
number of mixed vertices , and ch be \Vj(H, S) n (P(Vc) U Vc)\- By Lemma [6j the time 
complexity of computing Vj+i(S) from Vj(S) is 

^2 ((h H + m' H + c H )\og(n + c)). 
HeS,\H\=j 

By TheoremJTJ YlneS,\H\=j h H+m' H = 0(j(n-j)+jc) . It is also clear that Etfes,|H|=j c h = 
0(c). Therefore, the total time complexity is 0((j(n — j) + jc) log(n + c)). The correctness 
follows from the correctness proof of Lemma [6j □ 

Theorem 3 Vk(S) can be computed in 0(k 2 (n + c) log(n + c)) time. 

Proof. By Lemma [7J the total time complexity is Yli=i 0((i(n — i) + ic) log(n + c)) = 
0(fc 2 (n + c)log(n + c)). □ 
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4.2 Divide-and-Conquer Algorithm for Farthest-Site City Voronoi Diagram 

In this section we describe a divide-and-conquer approach to compute the farthest-site city 
Voronoi diagram J T V(S). Since there are n Voronoi regions in J r V(S) and each of them is 
associated with a site p G S, we denote such a region by TV(p, S). 

The idea behind this algorithm is as follows: To compute J-'V(S), divide S into two 
equally-sized sets Si and S2 = S\Si, compute ■FV(S'i) and J r V(S2), and then merge the two 
diagrams into J r V(S). Now, suppose we have already computed J r V(5i) and J r V(S < 2). Then, 
the edges of a Voronoi region FV(p, S) in FV(S) stem from three sources: i) contributed 
by J r V(5i), ii) contributed by J r V(<S , 2), and iii) contributed by two points, one in Si and 
the other in £2, that have the same distance to two farthest site. In fact, the union of all 
of the third kind of edges is Bq{Si, 1S2). Each connected component of Bc(Si, S2) is called 
a merge curve. A merge curve can be either a closed or open simple curve. 

If all the merge curves are computed, merging J r V(Si) and J-V{S2) takes time linear in 
the complexity of Bq{Si, S2). To compute the merge curves, we first need to find a point 
on each merge curve, and then trace out the merge curves from these discovered points. 

In order to compute a merge curve, we modify Cheong et al.'s divide-and-conquer al- 
gorithm [8] for farthest-polygon Voronoi diagrams in the Euclidean metric to satisfy our 
requirements. Given a set V of disjoint polygons, V = {Pi, P2, ■ ■ ■ , P m }, °f total com- 
plexity n, the farthest-polygon Voronoi diagram J r V('P) partitions the plane into Voronoi 
regions such that all points in a Voronoi region share the same farthest polygon in V . Let 
\P\ be the number of vertices of a polygon P G V and let \V\ be Ylp^v W\ = n - 

Their algorithm computes the medial-axis Ai(P) for each polygon P G V and refines 
FV{P 1 V) by Ai(P). Ai(P) partitions the plane into regions such that all points in a 
region share the same closest element of P, where an element is a vertex or an edge of 
P. In other words, for each point v G M 2 , M.{P) provides a shortest path between v and 
P. Therefore, the medial axes for J r V(P,V), with P G V, have the same function as the 
shortest path maps SVA4 P , with p G S in the city metric. By replacing V and Ai(P) with 
S and SVAip respectively, the divide-and-conquer algorithm of Cheong et al. [8] can be 
modified to compute J 7 V(S) with respect to the city metric. 

Cheong et al. [H] pointed out the bottleneck with respect to running time is to find 
for each closed merge curve a point that lies on it. In order to overcome the bottleneck, 
the authors use some specific point location data structures |9|12j . Let V be divided into 
two sets Vi and V2 = V \ V\, where \V\\ ~ l^l ~ \- Cheong et al. [Sj construct the 
point location data structures for TV(Vi) and J 7 V(T > 2)- For each polygon P G Vi and each 
vertex v G M.(P)nJ r V(P,Vi), they perform a point location query in FV{Vi) and J r V(T , 2) 
(likewise for each polygon P' G V%)- Each point location query requires O(logn) primitive 
operations, and each operation tests for 0(1) points and takes O(logn) time. Hence, one 
point location query takes 0(log 2 n) time. Since |-FV(7>i)| = |-FV(7>i)| = 0(n), merging 
TV{Vi) and TV(V 2 ) takes 0{n log 2 n) time. 



Since in our case |J r V(5i)| = |J 7 V(S , 2)| = 0(nc), we perform 0(nc) point location 
queries, each of which takes 0(log 2 nc) = 0(log 2 (n+c) 2 ) = 0(log 2 (n + c)) time. Therefore, 
merging J : V{S\) and FV{S2) takes 0{nc log 2 (n + c)) time. We conclude: 

Theorem 4 J r V(S') can be computed in 0(nc log n log 2 (n + c)) £ime. 

Proof. In the beginning, for each site p £ S, FV{{p}) is exactly SVM. V . Computing SVA4 P 
takes O(clogc) time [6], implying that computing J 7 V({p}), for all p £ S, takes 0(nc log c) 
time. Consider the merge process at some level i. The set S is divided into 2* subsets, and 
each of them contains at most n/2* sites. Therefore, the merging process at level i takes 
2 l ■ 0{n/2 l log 2 (n/2* + c)) = 0{nc log 2 (n + c)) time. Since there are logn levels, J 7 V(5) can 
be computed in O {nc logn log 2 (n + c)) time. □ 

5 Conclusion 

We contribute two major results for the /c th -order city Voronoi diagram. First, we prove 
that its structural complexity is 0(k(n — k) + he) and i?(n + kc). This is quite different 
from the 0(k(n — k)) bound in the Euclidean metric [llj. It is especially noteworthy that 
when k = n — 1, i.e., the farthest-site Voronoi diagram, its structural complexity in the 
Euclidean metric is 0(n), while in the city metric it is 0(nc). Secondly, we develop the first 
algorithms that compute the fc th -order city Voronoi diagram and the farthest-site Voronoi 
diagram. Our algorithms show that traditional techniques can be applied to the city metric. 
Furthermore, since the complexity of the first-order city Voronoi diagram is 0(n + c), one 
may think that the complexity the transportation network contributes to the complexity 
of the fc th -order city Voronoi diagram is independent of k. However, our results show that 
the impact of the transportation network increases with the value of k rather than being 
constant. 
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